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p£ : Abstract 



We numerically calculate the density profile and excitation spectrum of a two- 
species Bose-Einstein condensate for the parameters of recent experiments. 
We find that the ground state density profile of this system becomes unstable 
in certain parameter regimes, which leads to a phase transition to a new stable 
state. This state displays spontaneously broken cylindrical symmetry. This 
■ behavior is reflected in the excitation spectrum: as we approach the phase 

transition point, the lowest excitation frequency goes to zero, indicating the 
onset of instability in the density profile. Following the phase transition, this 
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I. INTRODUCTION 



Following the first observation of a two-species BEC by Myatt et al QXJ] there has been 
increased interest in the properties and applications of two-species Bose-Einstein conden- 
sates. Theoretically, two-species BECs are interesting because they allow the relative phase 
between the mean field wavefunctions to become manifest, whereas the phase of a single 
BEC cannot be observed. The original experiment of Myatt et al used the \F, rrif) = |1, — 1) 
and 1 2 , 2) hyperfine sublevels of 87 Rb in a cloverleaf trap. The two species formed slightly 
overlapping clouds which were offset relative to one another by gravity. Recent progress has 
also been made in the production of two-species BECs in which the two condensates share 
the same trap center, and are therefore less like two separate single condensates 0. This 
system uses the 1 1 , — 1 ) and \2, 1) sublevels of 87 Rb which have the same magnetic dipole 
moment to first order. 

In this paper, we consider the quasiparticle, or collective excitation, spectrum of this 
latter system We find an interesting dependence of the excitation frequencies which 

clearly shows the onset of instability in the density profile as particle number is increased, 
and a subsequent phase transition to a stable configuration. This behavior is due to the 
tendency of the two species to separate out like oil and water in certain parameter regimes |7| . 
As in the two-dimensional calculation of Ohberg et al || , our three-dimensional calculations 
show that this phase transition can take the form of a spontaneous breaking of cylindrical 
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symmetry. The onset of this instability as particle number is increased is indicated by an 
excitation frequency which goes to zero. Following a phase transition to a stable state, this 
frequency rises again to a constant value. 

Our work builds on previous studies of the excitations and stability properties of two- 
species BECs. Graham and Walls |3|] have analytically solved for the excitations of a binary 
phase two-species BEC in a trap, in the Thomas-Fermi limit and in the case where the 
parameters satisfy conditions leading to a pure binary phase two-species BEC. Busch et al 
Q have performed a variational calculation in which they compute the lowest lying excitation 
frequencies. Esry and Greene have numerically calculated the excitation spectrum of a 
two-species BEC in the in a TOP trap in which gravity separates the centers of the two 
condensates. They found that the repulsive interaction between the two clouds lead to 
simultaneous collective excitations of both species which were significantly different from 
the single condensate case. However, this system did not display the stability properties 
presented here. Pu and Bigelow [[],[| have solved for the case of mixed Rb-Na condensate 
under assumptions of cylindrical symmetry of the density profiles, and found that this system 
also displayed instabilities of the kind presented here, however, they have not investigated 
the effect of spontaneous breaking of cylindrical symmetry. 

II. GROUND STATE DENSITY PROFILE AND EXCITATION SPECTRUM 

For a two-species condensate, the second quantized grand canonical Hamiltonian is, in 
the position basis, 
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where ^ A (r) and ^s(r) are the field annihilation operators for species A and B respec- 
tively, fi A /B are the chemical potentials for the two species, m is the atomic mass (assumed 
here to be equal for species A and B), V A /b are the trapping potentials for the two species 
and the U AA /bb/ab are the scattering parameters (U pq = (A-Ka pq h 2 )/m, where cl aa / bb/ab are 
the scattering lengths between two atoms of species A, two atoms of species B and an atom 
of species A and an atom of species B). 

Setting ^ A = ijj + 5 A and = ip + 5b, where the ip = (^f\ are c-numbers and the <5's 

contain the operator dependence, we can expand Eq. (|l|) in the deviation operators 5 A /b- 
Assuming that the condensate is approximately in a coherent state allows us to discard the 
third and fourth order terms in the <5's, since these are small. The zeroth order term in the 5s 
gives the energy functional, whose stationary points subject to a particle number constraint 
are given by solutions to the Gross-Pitaevskii equation ||10|| . The first order term vanishes 
if the Gross-Pitaevskii equation is satisfied, and the second order term can be diagonalized 
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with a Bogoliubov transformation, giving the quasiparticle or excitation energies for the 
BEC. 

We use a variation of the well known basis set method [[R] to minimize the energy 



functional. This method consists in working in a basis of harmonic oscillator eigenfunctions 
for the bare trapping potential. We scale these basis vectors appropriately so that the 
expansion of the condensate due to repulsive atom-atom interactions is partially accounted 
for. The scaling factors are determined by a variational technique using Gaussian trial 
wavef unctions. This allows us to calculate using fewer basis vectors, or even, in certain 
geometries, to accurately scale out the dependence in one or two dimensions. 
We define a scaled basis as follows 



Ux, y, z) = (\ x \ y \ z ) 1/2 §i ( f , f , ±) (2) 

\ Ax Ay A z j 

where the $j(r) are normal modes of the free trap Hamiltonian and the A's are scaling 
factors determined by minimizing the energy functional with Gaussian trial wavef unctions. 
We have found that we can generally achieve reasonable accuracy with around two hundred 
basis vectors, although this depends on the geometry of the problem and the parameter 
regimes explored. 

In this basis, the energy functional becomes 
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where the oti and /3j's are the amplitudes in the various modes for species A and B respec- 
tively, H^/ B are the matrix elements of the non-interacting Hamiltonian p 2 /2m + V A /b{t) 
in the expanded basis and the V^^ AB ^ BB are the matrix elements of the two body potential: 

Vgb = J ^ 3 ^(r)^(r)^(r)^(r) (4) 

where p and q denote the species of atom (A or B). The stationary conditions of the energy 
functional subject to constraints of particle number conservation lead to the Gross-Pitaevskii 
equations for a two-species condensate |ll]||,|9j . 

We determine the ground state density profile by numerically minimizing the energy 
functional (|j) subject to constant particle number. We do this rather than the more usual 
and computationally simpler procedure of solving the Gross-Pitaevskii equation since any 
minimum of the energy functional is a solution to the Gross-Pitaevskii equation, but the set 
of solutions to the latter also include maxima and saddle points of the energy functional. In 
the particular case considered here, solutions to the Gross-Pitaevskii equation can become 
unstable as parameters are varied, and we need to take this into account. This instability 
is a result of the tendency of two-species condensates to form two separate clouds in certain 
parameter regimes. 

The excitation spectrum is then calculated by solving the Bogoliubov equations. In the 
discrete basis used here, these equations amount to solving the eigenproblem for a non- 
Hermitian matrix formed from the as and /3s. The two-species condensate case generalizes 
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readily from the single condensate case, for which details can be found in |]I2" |. The matrix 
to be diagonalized has the form: 



where 
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P BB and Q BB are defined as for P AA and Q AA respectively with A <->■ B and a <-> (5. The 
various symmetries which exist for given trapping geometries and parameter regimes can be 
used to simplify the problem by separating the matrix ([5]) into several smaller matrices. 



III. RESULTS 

We consider here equal trapping potentials for species A and B. This situation has 
been experimentally achieved for 87 Rb atoms in the \F, Mp) = |1, —1) and |2, 1) hyperfine 
sublevels which we label \A) and \B) JJ. We consider trap frequencies of f x = f y = 47/\/8 
Hz, f z = 47 Hz which are relevant to these experiments in a TOP trap. The low inelastic 
collision cross section for these species is consistent with scattering lengths claa/bb/ab which 
are all similar. Recent calculations give a A A = 5.68 nm and clbb = 5.36 nm [13]. We consider 
the three cases cl A b = 5.0 nm, 5.52 nm and 6.0 nm, which are within the range of expected 
values for this quantity, in order to show the onset of phase instability as a A B is increased. 

As shown by Ho and Shenoy 0, a binary condensate in the high particle number 
(Thomas-Fermi) limit can contain volumes in which only one species is present (giving 
two single particle phases, one for each species) and volumes in which both species coexist 
(binary phase) separated by phase boundaries. As an example, a condensate for which all 
parameters except the interspecies scattering length are the same for each species will exist 
as a single binary phase cloud only if U A b < U AA /bb- If we discard the Thomas- Fermi as- 
sumption, then this single species phase / binary phase picture is only approximate since the 
kinetic energy term in the Hamiltonian precludes the existence of sharp phase boundaries. 
However, a two-species BEC can still undergo phase transitions in the sense that a solution 
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to the Gross-Pitaevskii equation may become unstable as parameters are varied, and the 
condensate will then undergo collapse to some other stable solution. 

FIG. (0) shows a radial cross section through the density profile for clab — 5.0 nm and 
Na = N B = 15000 atoms. We see that in this case the two-species condensate exists as two 
highly overlapping clouds with no tendency to repel each other. The corresponding excita- 
tion spectrum is shown in FIG. (|| (b)), and FIG. (|| (a)) shows a single species condensate 
with a = y/dAAdBB- We see that, as predicted by Graham and Walls ||, the excitation 
spectrum of the two-species condensate undergoes a doublet splitting as compared to the 
case of a single species condensate. We have compared our results to the high atom number 
limit derived by Graham and Walls for a spherically symmetric two-species condensate and 
found agreement to within a few percent with less than two hundred basis vectors. 

FIG. (^|) shows what happens when clab is increased to ^/oaMIbb = 5.52 nm. We see 
that at high atom number, the lowest energy eigenvalue approaches closer to zero frequency 
than in the case of FIG. (Q), suggesting that we are near a region of phase instability in 
parameter space. 

Indeed, in the case clab — 6.0 nm, the single phase solution to the GP equation becomes 
unstable at a critical atom number of around 4000. The minimum of the energy functional 
is given by a different density profile, as shown in FIG. (f4|). It can clearly be seen that the 
cylindrical symmetry is spontaneously broken by the two condensates' mutual repulsion. 
This effect was also seen by Ohberg and Stenholm |§ in a two dimensional calculation. It 
illustrates the danger in assuming spherical or cylindrical symmetry when solving for the 
density profile of a binary condensate, since in this case such an assumption leads to a state 
which is unstable with respect to certain antisymmetric perturbations. 

The excitation spectrum for the parameters of FIG. (|4]) is shown in FIG. (|5|). We can 
clearly see that, for low atom numbers, the lowest (antisymmetric) excitation frequency goes 
to zero, suggesting the increasing instability of the condensate to antisymmetric perturba- 
tions. This is indeed the case, as is shown by the breaking of cylindrical symmetry seen in 
FIG. (Q). When N « 4000 is reached, the condensate undergoes a phase transition to an 
unsymmetric ground state and the lowest energy eigenvalue increases again to an asymptotic 
limit. 

One further interesting feature of this excitation spectrum is the persistence of excitations 
which look like the single condensate case even after the phase transition is reached - for 
example, the nearly horizontal excitation at 16.6 Hz which is hardly modified after the phase 
transition point. We interpret this as being due to the fact that, in the region above the 
phase transition the combined density profile of the two species looks very much like a single 
species condensate, even though the individual density profiles of each species are greatly 
modified (see FIG. (|j). The excitations in question are then interpreted as the normal single 
species type of excitations for this combined density profile. 

IV. CONCLUSION 

In certain parameter regimes, two-species Bose-Einstein condensates can show more 
complex density profile behavior than single condensates. In particular, solutions to the 
Gross-Pitaevskii equation can become unstable or metastable as parameters are adiabati- 
cally changed, leading to different ground state solutions. We have examined the way this 
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effect shows up in the excitation spectra of experimentally realistic two-species condensates 
and found that the onset of instability is heralded by an excitation frequency which tends 
to zero, as would be expected. The nature of the instability is consistent with the parity of 
the corresponding excitation. 
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FIGURES 

FIG. 1. Radial cross section through the density profile for the parameters clab = 5.0 nm 
and Na = Nb = 15000 atoms. The cross section is taken through the minimum of the trapping 
potential in the longitudinal direction. The mesh plot shows the maximum density for species 
A and B, and the line plots show the densities along the x and y axes of the combined density 
(upper plots) and species A and B densities (solid and dotted lines respectively). In this case, the 
scattering length clab is small enough compared to the two single species scattering lengths that 
the two species show no tendency to form separate clouds. 

FIG. 2. (a) Excitation spectrum for a single species BEC with clab = 5.52 nm. Parameters 
were chosen to approximate (b) as closely as possible, (b) Excitation spectrum for the parameters 
of FIG. (0). 

FIG. 3. Excitation spectrum for clab = 5.52 nm. 

FIG. 4. Density profile for clab = 6.0 nm. The top plot shows the densities for ./V = 3000 
and the bottom plot shows the densities for N = 4500. The spontaneous breaking of cylindrical 
symmetry is seen in the latter case. 

FIG. 5. Excitation spectrum for clab = 6.0 nm. The onset of phase instability can be seen 
at around N = 4000 and is characterized by the lowest excitation frequency approaching zero. 
Following a phase transition to a symmetry broken state, the lowest energy eigenvalue increases 
again, tending towards an asymptotic limit. 
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